To appear in The Astrophysical Journal 



Neutron-capture elements in the s- and r-process-rich stars: 
Constraints on neutron-capture nucleosynthesis processes 

Bo Zhang, Kun Ma, and Guide Zhou 

Department of Physics, Hebei Normal University, Shijiazhuang 050016, China; 

zhangbo@hebtu. edu. en 

ABSTRACT 

The chemical abundances of the very metal-poor double-enhanced stars are 
excellent information for setting new constraints on models of neutron-capture 
processes at low metallicity. These stars are known as s+r stars, since they show 
enhancements of both s-process and r-process elements. The observed abundance 
ratios for the double-enhanced stars can be explained by those of stars that were 
polluted by an AGB star and subsequently accreted very significant amounts of 
r-process material out of an AIC (accretion-induced collapse) or Type 1.5 super- 
nova. In this paper we present for the first time an attempt to fit the elemental 
abundances observed in the s- and r-rich, very metal-poor stars using a para- 
metric model and suggest a new concept of component coefficients to describe 
the contributions of the individual neutron-capture processes to double-enhanced 
stars. We find that the abundance ratios of these stars are best fitted by enrich- 
ments of s- and r-process material. The overlap factor in the AGB stars where 
the observed s-process elements were produced lies between 0.1 and 0.81. Tak- 
ing into account the dependence of the initial-final mass relations on metallicity, 
this wide range of values could possibly be explained by a wide range of core- 
mass values of AGB stars at low metallicity. The component coefficient of the 
r-process is strongly correlated with the component coefficient of the s-process 
for the double-enhanced stars. This is significant evidence that the r-process 
material in double-enhanced stars comes from an AIC or Type 1.5 supernova. 

Subject headings: nuclear reactions, nucleosynthesis, abundances— stars: AGB 
and post-AGB— supernovae: general 
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1. Introduction 

The elements heavier than those belonging to the iron peak are mostly made through 
neutron capture via two principal processes: the r-process (for rapid neutron-capture process) 
and the s-process (for slow neutron-capture process; Burbidge et al. 1957). During quiescent 
He burning, neutrons generated via the ^^C{a, ny^O (Cameron 1955; Reeves 1966; Mathews 
& Ward 1985) and ^^Ne(Q;, n)^^Mg (Cameron 1960) reactions are captured by heavy-element 
seed nuclei. The neutron densities are relatively low {Nn ~ 10^ cm~^), so nearly all possible 
/3-decays will have time to occur between successive neutron captures. Successful synthesis 
of heavier isotopes progresses along the "valley of /3 stability." This synthesis route is called 
the s-process, and it is responsible for the production of about half of the isotopes of the n- 
capture elements. The overabundances of elements heavier than iron observed at the surface 
of asymptotic giant branch (AGB) stars (Smith & Lambert 1990) clearly indicate that the 
s-process takes place during the AGB phase in the evolution of low- and intermediate-mass 
stars (0.8 < M/Mq < 8). The site or sites of the r-process are not known, although 
suggestions include the u-dnven wind of Type II supernovae (Woosley & Hoffman 1992; 
Woosley et al. 1994), the mergers of neutron stars (Lattimer & Schramm 1974; Rosswog et al. 
2000), AIC (accretion-induced collapse; Qian & Wasserburg 2003), and Type 1.5 supernovae 
(Zijlstra 2004). Extremely neutron rich nuclei are formed in a matter of seconds in the r- 
process, and it is not necessary to have preexisting heavy-element seed nuclei. The r-process 
is also responsible for about half of the abundances of solar system n-capture isotopes. Thus, 
n-capture elements can be composed of some pure r-process, pure s-process, and some mixed- 
parentage isotopes. As a result, when the solar system total abundances {tgs) are separated 
into contributions from the s-process {sgs) and the r-process (r^^), some elements are mostly 
contributed by the r-process, such as Eu, and some by the s-process, such as Ba. Therefore, 
Eu is commonly referred to as an "r-process element," and Ba, as an "s-process element." 

Most metal-poor stars in the Galaxy show that about 14% of those are carbon enhanced 
(hereafter CEMP stars; Cohen et al. 2005). These stars are often enriched in neutron-capture 
heavy elements. However, the enhancement of s-process and r-process elements varies from 
star to star. Observations of metal-poor stars with metallicities lower than [Fe/H]Ri —2.5 
enriched in neutron-capture elements have revealed the solar r-process pattern, while only 
two cases of highly r-process-enhanced stars (hereafter "r" stars), CS 22892-052 (Sneden 
et al. 1996, 2003) and CS 31082-001 (Cayrel et al. 2001; Hill et al. 2002), have been noted. 
Despite their considerable metal deficiency, these stars seem to have experienced an r-process 
that barely differs from the sum of r-processes that enriched the pre-solar nebula. This has 
led to suggestions that r-process production may be independent of the initial metallicity of 
the site, especially for the heavier n-capture elements (Z> 56; Cowan et al. 1995; Sneden et 
al. 1996; Norris et al. 1997; Sneden et al. 2000). 
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Although the material from which Population II stars form is not expected to contain 
significant s-process contributions, some stars including some subgiants are greatly enriched 
in carbon and s-elements (hereafter s stars; Norris et al. 1997; Hill et al. 2000). These are 
believed to be binary companions of initially more massive donor stars that have evolved 
through the thermally pulsing AGB phase and transferred material enriched in C and s- 
process elements onto the lower mass, longer lived secondary now observed. 

The discovery that several stars show enhancements of both r-process and s-process 
elements (Hill et al. 2000; Cohen et al. 2003) is puzzling, as they require pollution from 
both an AGB star and a supernova. Qian & Wasscrburg (2003) proposed a theory for the 
creation of "5 -|-r" -process stars. First some s-process material is accreted from an AGB 
star, which turns into a white dwarf. Later in the evolution of the system, the white dwarf 
accretes matter from the polluted star and suffers an AIC to a neutron star. The i/-driven 
wind produces an r-process, which also pollutes the companion. A possible problem these 
authors mention is the still uncertain nucleosynthesis in accretion induced collapse, which 
may or may not produce the r-process. Another possible s+r scenario is that the AGB 
star transfers s-rich matter to the observed star but does not suffer a large mass loss and at 
the end of the AGB phase; the degenerate core of the low-metallicity, high-mass AGB star 
may reach the Chandrasekhar mass, leading to a Type-1.5 supernova (Zijlstra 2004). Such 
supernovae can explain both the enhancement pattern and the metallicity dependence of the 
double-enhanced halo stars. 

There is another possibility for the double-enhanced scenario of the halo stars. In this 
picture, the formation of a binary system of low-mass stars was triggered by a supernova that 
polluted and clumped a nearby molecular cloud. Subsequently, the observed star, already 
strongly enhanced in r-process elements, receives large amounts of s-process elements from 
the initially more massive star that underwent the AGB phase, and turns into the double- 
enhanced star (hereafter "r+s" stars; Delaude et al. 2004; Barbuy et al. 2005; Gallino et al. 
2005; Ivans et al. 2005). Aoki et al. (2002) proposed that this scenario could possibly explain 
the formation of the double-enhanced stars, but the r-rich very metal-poor stars are so rare 
that it is implausible to explain the larger number of double-enhanced stars through even 
rarer ones. Previously, Barbuy et al. (2005) and Wanajo et al. (2006) suggested massive 
AGB stars (M Ki 8— I2M0) to be the origin of these double enhancements. Such a large- 
mass AGB star could possibly provide the observed enhancement of s-process elements in 
the first phase and explode or collapse providing the r-process elements, but the modeling 
of the evolution of such a large-mass metal-poor star is a difficult task, and the amount of 
the s-process material produced and its abundance distribution is still uncertain (Wanajo et 
al. 2006). 
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It is now widely accepted that the neutron exposure required to produce s-elements in 
AGB stars originates in some partial mixing of protons (PMP) from the envelope down into 
the C-rich layers resulting from the former intermittent operation of He burning during the 
thermal pulse phase of AGB stars. PMP activates the chain of reactions ^^C{p,'~fy^N{f3y^C{a,ny^0. 
The s-elements thus produced in the deep interior by successive neutron captures are sub- 
sequently brought to the surface by the third dredge-up (Gallino et al. 1998). Using the 
primary-like neutron source [^^C(q;, n)^^0] and starting with a very low initial metallic- 
ity, most iron seeds are converted into ^°^Pb. Thus, when the third dredge-up episodes 
mix the neutron-capture products into the envelope, the star appears s-enhanced and lead- 
rich. Therefore, if the standard PMP scenario holds, all s-process-enriched AGB stars with 
metallicities [Fe/H]<-1.3 are thus predicted to be Pb stars ([Pb/hs]>l, where hs denotes the 
'heavy' s-process elements such as Ba, La, and Ce), independent of their mass and metallicity 
(Goriely & Mowlavi 2000). 

The first three "lead" stars (HD 187861, HD 224959, and HD 196944) have been reported 
by Van Eck et al. (2001). At the same time, Aoki et al. (2001) found that the slightly 
more metal-deficient stars LP 625-44 and LP 706-7 are enriched in s-elements but cannot be 
considered as lead stars ([Pb/Ce]<0.4), in disagreement with the standard PMP predictions. 
A large spread of ^^C pocket efficiencies is proposed by Straniero et al. (2005) in order to 
explain the spreads of [Pb//is]. However, a physical explanation for the different -"^^C pocket 
strengths has not been yet found, so fundamental uncertainties currently exist in the models 
of AGB stars. 

The nucleosynthesis of neutron-capture elements in CEMP stars can be investigated by 
abundance studies of s-rich and r-rich stars. Recently, Cohen et al. (2003), by analysing 
the spectra of the double-rich star HE 2148-1247, concluded that the observed abundances 
could not be fitted by a scaled solar system s-process, r-process, and total solar inventory. A 
quantitative understanding of the origins of neutron-capture elements in the double-enhanced 
halo stars has so far been a challenging problem. Although some of the basic tools for this task 
were presented several years ago, the origins of the neutron-capture elements in the double- 
enhanced halo stars, especially r-process elements, are not clear, and the characteristics of 
the s-process nucleosynthesis in the AGB stars are not ascertained. Clearly, the study of 
elemental abundances in these objects is important for investigation of the origin of neutron- 
capture elements in the objects and in the Galaxy. 

In this paper we investigate the characteristics of the nucleosynthesis pathway that 
produces the abundance ratios of double-rich objects using the s-process parametric model 
and considering very significant amounts of enrichment of r-process material out of an AIC 
or Type 1.5 supernova(SN 1.5). The parametric model of double-enhanced stars is described 
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in Sect. 2. The calculated results are presented in Sect. 3, in which we also discuss the 
characteristics of the s- and r-process nucleosynthesis. Conclusions are given in Sect. 4. 

2. Parametric Model of Double-enhanced Stars 

There have been many theoretical studies of s-process nucleosynthesis in low-mass AGB 
stars. Unfortunately, however, the precise mechanism for chemical mixing of protons from the 
hydrogen-rich envelope into the ^^C-rich layer to form a ^^C pocket is still unknown (Busso 
et al. 2001). This makes it even harder to understand the particular abundance pattern 
of the s- and r-process elements found in carbon-rich, metal-poor stars. It is interesting to 
adopt the parametric model for metal-poor stars presented by Aoki et al. (2001), with many 
of the neutron-capture rates updated (Bao et al. 2000), to study which physical conditions 
reproduce the observed abundance patterns found in the stars. On one hand, this approach 
has the virtue of being model independent, in that it does not refer to detailed stellar 
evolution models. On the other hand, strictly speaking, our model does not represent exactly 
the situation proposed by current stellar s-process models. In the parametric model, in fact, 
we make use of an exponential distribution of neutron exposures via sequential irradiations 
(Howard et al. 1986). This representation would well describe a case in which the ^^C neutron 
source burns during thermal pulses. In the current stellar s-process model, instead, the ^^C 
neutron source is activated in radiative conditions in the periods between thermal pulses 
(Straniero et al. 1995), and its activation cannot be approximated by a simple exponential 
law of neutron irradiations (Gallino et al. 1998). There are four parameters in the parametric 
model on s-process nucleosynthesis. They are the neutron irradiation time At, the neutron 
number density the temperature Tg (in units of 10^ K) at the onset of the s-process, and 
the overlap factor r. Combining these quantities, we can obtain the neutron exposure per 
thermal pulse. At = NnVr^t, where Vt is the average thermal velocity of neutrons at Tg. 
The temperature and the neutron number density Nn are fixed respectively at reasonable 
values for the ^^C(Q;,n)^^0 reaction, Tg = 0.1 K and A'„=10^ cm^^, for these studies. Thus, 
the final s-process abundance distributions depend only on the neutron exposure At and 
overlap factor (Aoki et al. 2001). 

We explored the origin of the neutron-capture elements in the double-enhanced stars by 
comparing the observed abundances with predicted r- and s-process contributions. For this 
purpose, we propose that the ^th element abundance in the star can be calculated as 

Ni{Z) = CsNi, s + CrNi, .lOl^'^/^J, (1) 

where Z is the metallicity of the star, Ni^ s is the abundance of the ^th element produced 
by the s-process in the AGB star and ^ is the abundance of the ith element produced 
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by the r-process (per Si=10^ at Z—Zq), and Cg and Cr are the component coefficients that 
correspond to contributions from the s-process and the r-process, respectively. It should be 
noted that the s-process abundances in the envelopes of the stars could be expected to be 
lower than the abundance produced by the s-process in the AGB star, because the material 
is mixed with the envelopes of the primary (former AGB star) and secondary stars. 

Based on equation (1), we carry out the s-process nucleosynthesis calculation combined 
with the contribution of the r-process to fit the abundance profile observed in the double- 
enhanced stars, in order to look for the minimum in the four-parameter space formed by 
r, Ar, Cs, and C,.. The adopted initial abundances of seed nuclei lighter than the iron peak 
elements were taken to be the solar-system abundances, scaled to the value of [Fe/H] of the 
star. Because the neutron-capture-element component of the interstellar gas that formed 
very metal-deficient stars is expected to consist of mostly pure r-process elements, for the 
other heavier nuclei we use the solar system r-process abundances (Arlandini et al. 1999), 
normalized to the value of [Fe/H]. 

The ultra-metal-poor star CS 22892-052 merits special attention, because this star has an 
extremely large overabundance of n-capture elements relative to iron and very low metallicity 
with [Fe/H]= -3.1. Many studies (Cowan et al. 1999; Sneden et al. 1996, 1998, 2000; 
Norris et al. 1997; Pfeiffer et al. 1997) have suggested that the abundance patterns of the 
heavier (Z>56) stable neutron-capture elements in very metal-poor stars are consistent with 
the solar system r-process abundance distribution. This concordance breaks down for the 
lighter neutron-capture elements in the range of 40<Z<56 (Sneden et al. 2000). Zhang et al. 
(2002) have reported that when the abundances of the lighter elements in CS 22892-052 are 
multiplied by a factor of 1/0.4266, the abundance distributions obtained for both heavier and 
lighter neutron-capture elements are in accordance with the solar system r-process pattern. 
This star could well have abundances that reflect the nucleosynthesis of a single supernova 
(Fields et al. 2002), so the adopted abundances of nuclei N^^ ^ in equation (1) are taken to 
be the solar system r-process abundances (Arlandini et al. 1999) for the elements heavier 
than Ba; for the other lighter nuclei we use solar system r-process abundances multiplied by 
a factor of 0.4266. 



3. Results and Discussion 

Using the observed data in 12 sample stars (Barbuy et al. 2005; Hill et al. 2000; Cohen 
et al. 2003; Aoki et al. 2001, 2002; Ivans et al. 2005), the model parameters can be obtained. 
The results of the neutron exposures, overlap factor, and the component coefficients are 
listed in table 1. 
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Figure 1 shows our calculated best-fit results. For most stars , the curves produced 
by the model are consistent with the observed abundances within the error limits. The 
agreement of the model results with the observations provides strong support to the validity 
of the parametric model adopted in this work. 

In the AGB model, the overlap factor r is a fundamental parameter. The overlap factor 
deduced for the double-enhanced stars lies between 0.1 and 0.81. Aoki et al. (2001) have 
reported an overlap factor of r~0.1 for metal-deficient stars LP 625-44 and LP 706-7; our 
calculated results of the two stars are close to their values. For the third drcdgc-up and the 
AGB model, several important properties depend primarily on the core mass (Iben 1977; 
Groenewegen & de Jong 1993; Karakas 2002). In the core-mass range 0.6 < Mc < 1.36, 
an analytical formula for the AGB stars was given by Iben (1977) showing that the overlap 
factor increases with decreasing core mass. Gombing the formula and the initial-final mass 
relations (Zijlstra 2004), we can obtain the overlap factor as a function of the initial mass 
and metallicity. In AGB stars with initial mass in the range M = 1.0 — 4M0, the core mass 
Mc lies between 0.6 and 1.4M0 at [Fe/H]~-2.5. According to the formula by Iben (1977), the 
corresponding values of r would range between 0.8 and 0.16. In an evolution model of AGB 
stars, a small r may be realized if the third dredge-up is deep enough for the s-processed 
material to be diluted by extensive admixture of unprocessed material. In fact, more recent 
AGB models showing the third dredge-up for low-mass stars at solar metaUicity have found 
(Gallino et al. 1998), for core masses around 0.7 M0, much lower overlap factors, down to 
r ~ 0.4, than the value of 0.8 given by Iben (1977). Karakas (2002) and Herwig (2000, 2004) 
have found that the third dredge-up is more efficient for the AGB stars with larger core 
masses, confirming the low values of r obtained by Iben (1977) in these cases. Taking into 
account the core-mass dependence, the wide range of r-values for the double-enhanced stars 
that we obtain could possibly be explained by a wide range of core-mass values of AGB stars 
at low metallicity. 

The neutron exposure per pulse Ar is also a fundamental parameter in the AGB model. 
The neutron exposure per pulse Ar deduced for the double-enhanced stars lies between 0.45 
and 0.88 mbarn"^. Aoki et al. (2001) have reported a neutron exposure per pulse Ar ~ 
0.71 mbarn"^ for the metal-deficient star LP 625-44 and Ar ~ 0.80 mbarn"^ for LP 706-7; 
our calculated results for these two stars are close to their values. In the case of multiple 
subsequent exposures the mean neutron exposure is given by tq = — Ar/lnr. For the s- 
process nucleosynthesis in the AGB stars, the final abundance distributions depend mainly 
on the mean neutron exposure tq. Basically, for a given neutron exposure one can obtain a 
higher mean neutron exposure simply by increasing the overlap factor, because there is more 
material subsequently exposed to neutron fluxes, which favors the production of the heavier 
elements, such as Pb. This implies that the formation of a Pb star ([Pb/Ba]>1.0) could be 
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explained by a transfer of matter from the low-mass AGB star (M<3.OM0). However, one 
should remember that actually in the model with the ^^C pocket burning radiatively it is 
mostly the amount of ^^C that determines [Pb/Ba]. Introducing a large range of the ^"''C 
amount at a given metallicity also proved to be effective for explaining the large scatter of 
[Pb/hs] observed in the metal-poor AGB stars (Straniero et al. 2005). 

The presence of s-process elements, along with large enhancements of carbon, suggests 
that a mass-transfer episode from a former AGB star in a binary system took place. Thus, 
one major goal is to find an astrophysical scenario, associated with an AGB star in a binary 
system, in which the r-proccss might also occur. For a binary system with small orbital 
separation, the matter accreted from the initial massive star is large (Boflin & Jorissen 1988), 
so the level of both s-enrichment and r-enrichment of the observed star should be high if the 
r-process also occurs in the binary system (e.g. AIC or SN 1.5). We compare the strength 
of the r-process elements parameterized by with the strength of the s-process elements 
parameterized by Cs for the double-enhanced stars. A surprising result shown in Figure 2 
is the strong correlation between and Cg. This implies an increase of s-process matter 
accreted from the AGB star with increasing r-process matter accreted from the AIC or SN 
1.5. The correlation obtained from this work is significant evidence for the s+r scenario. 

Based on the dependence of the initial-final mass relations on metallicity, Zijlstra (2004) 
speculates that the metal-poor AGB stars with initial mass in excess of 3 — 4Mq may cause 
the degenerate cores to reach the Chandrasckhar mass, leading to type-1.5 supernovae; this 
process would explain the enhancement pattern of the double-enhanced halo stars produced 
by small overlap factors (i.e. r~ 0.1 — 0.2). Since AIC is thought to occur only under 
somewhat restricted conditions, i.e., both the white dwarf mass and the mass transfer rate 
must be high enough, Bailyn & Grindlay (1990) assumed that the initial masses leading 
to heavy white dwarfs required for AIC would be in the range of 4:Mq < M < SMq and 
speculated that the rate of the occurrence of the AIC process in the Galaxy is no more 
than ~10~^yr~^. This rarity seems to be in confiict with the substantial fraction of double- 
enhanced stars (~30%) among all the CEMP— s stars currently observed. Cohen et al. 
(2003) proposed that AIC may be caused in binary systems with moderately massive AGB 
stars (M>3M0) and small initial orbital separation. Because the core mass of the AGB star 
is remarkably large at low metallicity (Zijlstra 2004), we propose that the binary systems 
with lower mass AGB stars (M<3— 4M0) may also cause AIC, which can increase the rate 
of the occurrence of AIC and explain the enhancement pattern of the double-enhanced halo 
stars produced by large overlap factors (r ~ 0.3 — 0.8). In this case, the orbital separation 
must be small enough to allow for capture of a sufficient amount of material to create the 
formation of AIC. The initial-final mass relation fiattens at higher metallicity (Zijlstra 2004), 
and the degenerate cores of high-metallicity AGB stars are smaller than those of the low- 
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metallicity stars, so the formation of AIC or SNe 1.5 is more difficult in fiigfi-metallicity 
binary systems, wfiicli could explain the upper limit of the metallicity ([Fe/H]<-2.0) for the 
observed double-enhanced stars. 

The binary systems with low-mass AGB stars (M<3Mo) and large initial orbital sep- 
aration could not cause AIC, because the white dwarf accretes matter insufficiently from 
the polluted star and the companion is polluted only by the former AGB star, which can 
explain the formation of s stars. For example, the r-process component coefficients of HD 
196944 and CS 30301-015 are 0.6 and 0.9, respectively, which implies that the two stars are 
s stars. The s-process component coefficients of these stars are smaller than those of the 
other 10 stars. This fact illustrates that there is less accreted material of these stars from 
the former AGB stars, and they should be formed in the binary systems with larger initial 
orbital separations. In addition, the s-element patterns of these two stars are matched by 
models with r=0.44 and 0.34, which are large enough to be consistent with the idea that the 
AGB star polluting these s stars was of low mass. 

The r star could possibly be formed in a close binary system, favoring an early mass 
transfer from the massive companion in its RGB(red giant branch) phase. In this case, the 
massive companion should produce no 5-elements, and its final explosion produces r-elements 
(Barbuy et al. 2005). The r star could be also formed in a molecular cloud that was polluted 
by a single supernova or a few supernovae (Shigeyama & Tsujimoto 1998; Fields et al. 2002). 
The halo star with perhaps the most striking r-process composition is CS 22982-052; its 
r-process abundance ratios are 

[^1.1.7^^.50(^)e. (2) 

This indicates Cr~50, in the range of component coefficients C^ obtained for the double- 
enhanced stars. 

As an example, it is interesting to investigate a possible explanation of the parameters 
obtained for a double-enhanced star using the s-|-r scenario. The s-process abundances of 
HE 2148-1247 are a result of pollution from the dredged-up material of the former AGB 
star. The measured [s/Fe] refers to the average s-processed material of the AGB star after 
dilution by mixing with the envelope of the presumed unevolved companion that is now the 
extrinsic star. The parameters deduced for HE 2148-1247 are r=0.1, At=0.88 mbarn"^, 
and Cs — 0.0045. Adopting the analj^ical formula given by Iben (1977) and the initial- 
final mass relations (Zijlstra 2004), we can know that the initial mass of the former AGB 
star is about 3—4 M©, which hes in the range of M ~ 3 — I2M0 reported by Cohen et al. 
(2003). At some time on the AGB, the convective He-shell and the envelope of the giant will 
be overabundant in heavy elements by factors f shell and fenv, 1, respectively, with respect 
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to solar system abundances normalized to the value of [Fe/H]. The approximate relation 
between /e„^,, i and f shell is 

J env, 1 ~ J shell: V'-'/ 

where AMdr is the total mass dredged up from the He shell into the envelope of the AGB 
star and is the envelope mass of the AGB star. For a given s-process element, the 
overabundance factor fenv, 2 in the future s+r star envelope can be approximately related 
to the overabundance factor fenv, 1 by 

_AM2 _ AM2 AMdr , . 

J env, 2 ~ J env, 1 ~ j^j^ j^^p J shell: [^) 

where AM2 is the amount of matter accreted by the future s+r star and M| is the envelope 
mass of the star. The component coefficient Cg is computed from the relation 

^ ^ fenv, 2 _ AM2 AMar . . 

' f shell ~ M| Mf ■ ^ ^ 

The value of f shell of course depends on the heavy element considered. For Pb, the factor 
f shell of HE 2148-1247 derived in the present study is about 1.6x10^, which lies in the 
6.3 X 10^ < f shell < 2.5 X 10^ range for the metal-poor AGB stars with metallicities from 
[Fe/H]=-2.0 down to -3.0 calculated by Busso et al. (2001). For the sake of simplicity, 
no continuous variation of the factor during the AGB phase has been considered. Taking 
AMdr = OMMq (Busso et al. 1992), AM2 ^ -0.2Mi (Steinhardt & Sasselov 2005), and 
M| «i O.5M0 (Boffin & Jorissen 1988), the value of Cs — 0.004 is deduced from equation 
(5), which is close to the value for HE 2148-1247. Because of the uncertainties of mass-loss 
rates and our poor knowledge of how and when mass transfer phenomena occur, we do not 
claim that this is the only or even the best understanding of the parameters. 

We discuss the uncertainty of the parameters using the method presented by Aoki et 
al. (2001). Figure 3 and 4 show the calculated ratios log(Pb/Ba) and log(Ba/ls), where 
Is denotes the 'fight' s-process elements (Sr and Y), as a function of the neutron exposure 
At in a model with r=0.1 and a function of overlap factor r with a fixed neutron exposure 
Ar=0.88 mbarn"^. These are compared with the observed ratios from HE 2148-1247. There 
is only a region of overlap in Figure 3, At =0.88±0.08 (1 a) mbarn"^, in which both the 
observed ratios log(Pb/Ba) and log(Ba/ls) can be accounted for. The bottom panel in 
Figure 3 displays the reduced value calculated in our model and there is a minimum, 
with — 2.36, at At — O.SSmbarn~^ with a 1 cr error bar ±0.08 mbarn"^; the neutron 
exposure is constrained quite well. The ratios log(Pb/Ba) and log(Ba/ls) are insensitive to 
the overlap factor r and allow for a wider range, < r < 0.28. The uncertainties of the 
parameters for the star HE 2148-1247 are similar to those for metal-poor stars LP 625-44 
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and LP 706-7 obtained by Aoki et al. (2001). For most sample stars, the uncertainties of 
the parameters are smaller than those for HE 2148-1247. 

Our model is based on the observed abundances of the double-enriched stars and nucle- 
osynthesis calculations, so the uncertainties of those observations and measurement of the 
neutron-capture cross sections will be involved in the model calculations. We note from table 
1 that for three stars (HE 2148-1247, LP 625-44, and CS 29497-030), the reduced are 
about 2. The probability that could be this large as a result of random errors in the mea- 
surement of the neutron-capture cross sections and the abundances of the neutron-capture 
elements is less than 2%. We find that all these uncertainties cannot explain the larger er- 
rors of neutron-capture elements, such as Zr in HE 2148-1247 and Y in CS 29497-030. This 
implies that the understanding of the true nature of the s-process or r-process is incomplete 
for at least some of these elements (Travaglio et al. 2004). 

4. Conclusions 

A main result of this work is the wide range of r-valucs for the double-enhanced stars. 
This can be explained when varying the initial mass of the AGB stars. We also find a 
strong correlation between the strength of the r-process elements and the strength of the 
s-process elements for the double-enhanced stars. This correlation is significant evidence 
for the scenario in which both the r- and s-enhancements of s-\-r stars are produced in 
the binary system. Since the formation of the AIC or SN 1.5 is more difficult for a binary 
system with high metallicity, the upper limit of the metallicity ([Fe/H]<-2.0) for the observed 
double-enhanced stars could be explained. 

Based on the discussion above, the double-enhanced stars belong to binaries. If these 

stars are found to be nonbinaries, they may also have been binaries initially, ceasing to be 
binaries due to the explosion of the massive companion. Overall, the number of CEMP stars 
is still small, especially for the CEMP-r stars, making proposed enrichment scenarios difficult 
to explicitly test. To underpin these studies, accurate abundance analyses for similar s- and 
r-process-rich, metal-poor, carbon-enhanced stars are required. Obviously, a more precise 
overlap factor-core mass law and initial-final mass relations from detailed models of AGB 
stars are needed. Further abundance studies of neutron-capture-rich stars will reveal the 
characteristics of the s- and r-processes at low metallicity, such as their metallicity and mass 
dependence, and the history of enrichment of neutron-capture elements in the early Galaxy. 

We are grateful to the referee for very valuable comments and suggestions that improved 
this paper greatly. This work has been supported by the National Natural Science Foundation 
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Fig. 1. — Best fit to observational results of metal-deficient stars. The black circles with 
appropriate error bars denote the observed element abundances, and the solid lines represent 
predictions from s-process calculations considering r-process contribution. 
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Fig. 2. — Correlation between the s-process component coefficients and the r-process 
component coefficients for the double-enhanced stars. The fit to the relation is shown as 
a solid line. 
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Fig. 3. — Abundance ratios log (Pb/Ba) (top), log (Ba/ls) (middle), and reduced (bot- 
tom), as a function of the neutron exposure per pulse, At, in a model with overlap factor 
r — 0.1. Solid curves refer to the theoretical results, and dashed horizontal lines refer to the 
observational results with errors expressed by dotted lines. The shaded area illustrates the 
allowed region for the theoretical model. 
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Fig. 4. — Same panels as those in Fig. 3, but as a function of the overlap factor r. 
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Table 1. Observed abundance ratios and the derived parameters for s- and r-rich stars 



Star 


[Fe/H] 


[Ba/Fe] 


[Eu/Fe] 


[Pb/Ba] 




r 


At 
(mbarn"^) 


To 

(inbarn"^) 










CS 22948-27 


-2.47 


2.26 


1.88 


0.46 


0, 


.37 


0.61 


0.614 


0.0033 


65.9 


0, 


.402859 


CS 29497-34 


-2.90 


2.03 


1.80 


0.92 


0, 


.61 


0.53 


1.072 


0.0034 


57.3 


1, 


.100821 


HE 2148-1247 


-2.30 


2.36 


1.98 


0.76 


0, 


.10 


0.88 


0.382 


0.0045 


67.4 


2, 


.360219 


CS 29526-110 


-2.38 


2.11 


1.73 


1.19 


0, 


.79 


0.64 


2.715 


0.0040 


50.6 


0, 


.580903 


CS 22898-027 


-2.25 


2.23 


1.88 


0.61 


0, 


.42 


0.77 


0.888 


0.0035 


67.9 


1, 


,225389 


CS 31062-012 


-2.55 


1.98 


1.62 


0.42 


0, 


.32 


0.71 


0.623 


0.0018 


37.3 


0, 


.883710 


CS 31062-050 


-2.32 


2.30 


1.84 


0.60 


0, 


.45 


0.71 


0.889 


0.0039 


60.6 


1, 


.027867 


HD 196944 


-2.25 


1.10 


0.17 


0.80 


0, 


.44 


0.45 


0.548 


0.0006 


0.6 


0, 


,586861 


CS 30301-015 


-2.64 


1.45 


0.2 


0.25 


0, 


.34 


0.54 


0.501 


0.0005 


0.9 


1, 


,118641 


LP 625-44 


-2.71 


2.74 


1.97 


-0.19 


0, 


,16 


0.69 


0.377 


0.0045 


76.2 


2, 


,111419 


LP 706-7 


-2.74 


2.01 


1.40 


0.27 


0, 


,10 


0.82 


0.356 


0.0017 


17.5 


0, 


,846228 


CS 29497-030 


-2.57 


2.32 


1.99 


1.33 


0, 


,81 


0.61 


2.895 


0.0060 


86.4 


2, 


,441984 



